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Abstract 

We propose a mean field theory for interfaces growing according 
to the Kardar-Parisi-Zhang (KPZ) equation in 1+1 dimensions. The 
mean field equations are formulated in terms of densities at different 
heights, taking surface tension and the influence of the nonlinear term 
in the KPZ equation into account. Although spatial correlations are 
neglected, the mean field equations still reflect the spatial dimension- 
ality of the system. In the special case of Edwards- Wilkinson growth, 
our mean field theory correctly reproduces all features. In the presence 
of a nonlinear term one observes a crossover to a KPZ-like behavior 
with the correct dynamical exponent z = 3/2. In particular we com- 
pute the skewed interface profile during roughening, and we study the 
influence of a co-moving reflecting wall, which has been discussed re- 
cently in the context of nonequilibrium wetting and synchronization 
transitions. Also here the mean field approximation reproduces all 
qualitative features of the full KPZ equation, although with different 
values of the surface exponents. 



1 Introduction 



Since many years the physical properties of deposition-evaporation processes 
on a planar surface have been studied theoretically by analyzing appropriate 



stochastic growth models that capture the essential features of the experi- 
mental realm [1]. In most of these models the configuration of the growing 
surface is described by a height variable h(x, t) that yields the height of the 
interface between deposited layer and gas phase above point x of the sub- 
strate at time t. Starting with a certain initial configuration, the interface 
then evolves according to certain stochastic rules. 

Depending on the specific dynamic rules for deposition and evaporation 
and their symmetries, the temporal evolution of the interface may be de- 
scribed on a coarse-grained scale by a stochastic differential equation, one of 
the simplest and more general one being the celebrated Kardar-Parisi-Zhang 
(KPZ) equation [2] 

^h(x, t) = vq + DV 2 h(x, t) + ^ [Vh(x, t)} 2 + £(x, t). (1) 

Here vo is the average growth velocity which can be set to zero in a co-moving 
frame, the Laplacian accounts for surface tension of the interface, and t) 
is an uncorrelated white Gaussian noise generated by the stochastic nature 
of deposition and evaporation. Moreover, the nonlinear term (V/i) 2 it's the 
simplest one which breaks the invariance under reflections h — > —h. 

As many models for interface growth, the KPZ equation exhibits dy- 
namic scaling, i.e., starting with a flat configuration the interface width 
w(t) = \J (h 2 ) — (h) 2 (where (•) denotes average over space and ensemble 
realizations) first increases as a power law w(t) ~ i 7 until it saturates in 
a finite system of linear size L at a stationary value lOgtat ~ L a . The 
crossover from a roughening to a stationary state is described by the well- 
known Family- Vicsek scaling form [3] 

w(t) ~ Fg(t/L*), (2) 

where g is a universal scaling function and z = a/7 is the dynamical expo- 
nent. 

The width is actually related to the second moment of the height dis- 
tribution profile Pi(h,t), which is defined as the normalized probability to 
find the interface at a randomly chosen lattice site at height h. Clearly, 
the height distribution contains much more information about the interface 
morphology than the width alone. As for the width, dynamic scaling implies 
a scaling form for the height distribution which in a co-moving frame may 
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be written as 

P L (h,t) = t~~<f(h/t\t/L z ) 



(3) 



Obviously, both the critical exponents a and 7 and the shape of height 
profile during roughening or after saturation reflect the symmetries of the 
growth process under consideration. The simple case of invariance under 
the reflection h — > —h can be studied by imposing A = in Eq. lfl|). In this 
case the linear Edwards- Wilkinson (EW) equation [4] is recovered, and the 
critical growth exponents in 1+1 dimensions take the values a. = 1/2 and 7 = 
1/4. Moreover, the height profile of 1+1-dimensional EW processes is known 
to be a simple Gaussian distribution both in the dynamically roughening 
phase as well as in the stationary state. 

In more realistic growth models, where nearest neighbour interactions 
play a role in the dynamics of the growing interface, reflection symmetry is 
broken and the nonlinear KPZ term has to be taken into account. In 1+1 
dimensions such a term is known to be relevant in the renormalization group 
sense. Therefore, even when the reflection symmetry is weakly violated 
(i.e., if A is small), the scaling behavior of an infinite system will eventually 
cross over from EW to KPZ scaling, the latter being characterized by the 
exponents a = 1/2, 7 = 1/3, and z = a/7 = 3/2. 

With a non-symmetric term being present there is no longer any reason 
for the height distribution Pi(h,t) to be symmetric with respect to h. Al- 
though in 1+1 dimensions a KPZ interface of a finite system after saturation 
still happens to be symmetric and Gaussian (see e.g. [5]), the profile of a 
roughening KPZ interface before saturation is indeed skewed [6], reflecting 
the asymmetry of the nonlinear term. In what follows we therefore restrict 
ourselves to the roughening process before saturation, i.e., t <C L z , regard- 
ing a virtually infinite system at finite times. Formally speaking, this can 
be achieved by taking the thermodynamic limit L — > 00 before the time- 
asymptotic limit t — > 00 is carried out. In particular, we are interested in 
the scaling function f(z) := f(z, 0), which renders the rescaled shape of the 
skewed profile after sufficiently long time (see left panel of Fig. U). 

The function f(z) is known to be universal, i.e., the asymptotic shape 
of the skewed profile is fully determined by the underlying KPZ field theory 
and does not depend on the microscopic details of the model. It has been 
suggested that the finite-time height distribution, especially the form of its 

lr The scaling functions g and / are related by g 2 (u) = J °° f(u,v)v 2 dv — 

[f™ f( v > u ) vdv ] 2 - 
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Figure 1: Skewed interface profile. Left panel: Numerically determined profile 
f(z) :— f(z,0) of a roughening KPZ interface in 1+1 dimensions for A < (solid 
line) compared to a Gaussian distribution of the same width (dashed line). Right 
panel: Effective exponents rj e ff(z) (see text) in comparison with the numerical 
estimates r] + = 1.6(2) and = 2.4(2) reported in Ref. [10] (representing the 
error bars as dashed boxes) and the theoretical predictions for the PNG model [8] 
(marked by full horizontal lines). 

tails, is approximately given by a stretched exponential 

P L (h,t)ocexp[-fi(\h-{h)\/t^] i«L 2 , (4) 

meaning that f(z) ~ exp(—[i\z\ v± ). Here [x is a metric factor while the expo- 
nents r\± refer to the two different tails of the distribution with ±X(h— (h)) > 
0. Because of the skewness both exponents are expected to be different. An 
argument based on a replica scaling analysis of directed polymers [7] , whose 
free energy fluctuations corresponds to the height fluctuations of a KPZ 
interface, suggests the value rj + = 3/2. 

As a breakthrough, Prahofer and Spohn have shown recently [8] that the 
finite-time rescaled height profile of the polynuclear growth model (PNG) 
[9], a model which is believed to belong to the KPZ universality class, equals 
the Gaussian orthogonal ensemble (GOE) Tracy- Widom distribution. This 
immediatly leads to 

/7+ = 3/2 , r?_=3. (5) 

Numerical simulations reported in literature concerning both directed poly- 
mers and KPZ lattice models [10] give r/ + = 1.6(2) and = 2.4(2), the 
latter value being not in agreement with theoretical predictions. However, 
since these results were obtained more than a decade ago the numerical 
precision was limited. Performing similar simulations using the so-called 
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Figure 2: Snapshots of a roughening interface in a reference frame where the 
asymptotic velocity is zero. The figure shows a free interface with A<0 (upper 
panel) compared to interfaces confined by a co-moving lower wall in the two cases 
A<0 (middle panel) and A>0 (lower panel). Simulations have been performed using 
the single step model (see Sec. 2) and snapshots have been taken after 2048 time 
steps. 



single step model [11] (see Sec. 2) we measured the effective exponent 
Vcffiz) = \ n f(z) lz ^ n f( z ) which according to Eq. (|1J) should converge to r/± 
as \z\ — > oo. As it can be seen in the right panel of Fig. ^ our numerical re- 
sults are in agreement with theoretical predictions for the PNG model while 
being incompatible with the previous estimate for of Ref. [10]. Clearly, 
the center of the scaling function (i.e. small values of z) is described only 
approximately by Eq. (|4"|). 

Looking at a snapshot of a roughening interface in 1+1 dimensions, it 
is almost impossible to recognize the influence of the KPZ nonlinearity by 
naked eye. Its influence, however, is much more pronounced in the presence 
of a hard-core wall. The wall is fixed in a frame where the asymptotic 
velocity of the interface vanishes and interacts with the interface solely by 
preventing excursions to negative heights. As can be seen in Fig. in 
presence of such a wall one can easily appreciate the dramatic difference 
emerging when the sign of A is changed 2 . Surprisingly, for A < the interface 
touches the wall only occasionally, while a high density of contact points is 

2 Alternatively one may compare a lower and an upper wall while keeping the sign of A 
fixed. 
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observed for A > 0. 

The properties of a KPZ interface close to a reflecting wall has been 
studied recently in the context of nonequilibrium wetting [12], where the 
interface describes a wetting layer on a planar substrate. Upon varying the 
average growth velocity vo the interface undergoes a depinning transition 
between a pinned phase, in which portions of the interface remains attached 
to the wall, to a depinned phase, where the interfaces detaches entirely and 
starts moving upwards. At the critical point, where the asymptotic interface 
velocity is zero, a second order phase transition takes place and various 
scaling laws can be singled out. In addition, the case A < describes the 
critical properties of most synchronization transitions in spatially extended 
chaotic systems [13,14]. 

Previous numerical simulations suggested that the temporal decay of the 
density po(t) of contact points, where the interface touches the wall, obeys 
the power law [15, 16] 

r 1.1(1) if A < 

p (t)~ t - e , e^l 3/4 if A = , (6) 
[ 0.22(2) if A > 

where the exponent 3/4 can be obtained from a transfer matrix calculation 
[12]. Moreover, a hyperscaling relation observed in simulations starting with 
a single pinned site [17] suggests the rational value 6 = 7/6. Obviously, the 
different values of the exponents reflect the asymmetry of the nonlinear term 
with respect to reflections h —* —h. Moreover, the pronounced numerical 
variation of 6 by a factor of 5 explains why the snapshots in Fig. [2] are so 
strikingly different. 

Interestingly, the profile of a roughening KPZ interface next to a wall 
cannot be described in terms of an appropriate generalized GOE Tracy- 
Widom distribution because of emerging nonlinear term. More generally, 
the critical behavior of such a bounded growth process in 1+1 dimensions 
is not easily accessible by analytical means. For example, renormalization 
group techniques fail either due to the presence of a strong-coupling fixed 
point unaccessible by perturbative approaches in the case A < [18] or due 
to essential singularities arising for A > 0. Therefore, the primary aim of 
the present paper is to discuss this case within a suitable mean field approx- 
imation. The mean field theory to be constructed should incorporate the 
asymmetry caused by the nonlinear term and should render a skewed height 
distribution with similar properties as in the full model. Although the mean 
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Figure 3: Single step model in 1+1 dimension. Left panel: The simulation usually 
starts with a flat interface, realized as a horizontal sawtooth pattern. On selecting 
a site with a local minimum a diamond (rhombus) is deposited with probability p, 
flipping up the interface by two units. Similarly, if the selected site happens to be 
a local maximum, a diamond may evaporate with probability 1 — p, flipping the 
interface downward by two units. Right panel: For p > i the interface roughens 
and propagates upwards. 



field theory ignores space, it should not resemble a naive infinite-dimensional 
limit (where a KPZ interface is always smooth), instead it should reflect to 
some extent the dimensionality of space in the thermodynamic limit L — > oo. 
Moreover, the desired theory should be as simple as possible and exactly 
solvable. In the following sections we propose and solve a mean field theory 
which meets these requirements. 



2 Mean field equations 

The mean field equations proposed here are inspired by a particular model, 
the so-called single step model (SSM) [11], which is probably the simplest 
and most compelling lattice model for KPZ-type interface growth. 

In the single step model the growing interface is represented by a set of 
integer heights fi; £ N residing at the sites i = 1 . . . L of a one-dimensional 
lattice of length L with periodic boundary conditions, obeying the restriction 

n i+1 -m = ±l. (7) 

The interface evolves in time according random sequential updates as fol- 
lows: At each sub-time step dt = 1/L a site i is chosen at random. If the 
interface has a local minimum at site i (i.e., < ni±\) the height is 
increased by 2 with probability p £ [0, 1] . This update can be pictured as 
depositing a diamond (see Fig. EJ, transforming a local minimum into a lo- 
cal maximum. Similarly, if the selected site happens to be a local maximum 
(rij > rii±i) the height rtj is decreased by two units with probability 1 — p. 
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For p = 1/2 the propagation velocity of the interface is zero and the 
evolution rules satisfy detailed balance, as described by the EW equation. 
For p 7^ 1/2, however, the propagation velocity is nonzero, depending on 
the roughness and the average slope of the interface. In this case the SSM 
exhibits KPZ growth with A being equal to \ — p. 

By identifying upward segments nj+i — n; L = 1 with particles and down- 
ward segments n^+i — n% = — 1 with vacancies, the single step model can be 
mapped exactly onto a partially asymmetric exclusion process (ASEP) [19] 
of diffusing particles with density 1/2. Since the ASEP is known to evolve 
towards an uncorrelated product state with a current j = p/2 — 1/4, it is 
immediately clear that the propagation velocity of the interface tends to 

= lim v(t) = p - I- (8) 

t^oo z 

as t — > cx). 3 Moreover, the mapping to the ASEP allows one to solve the 
model via Bethe ansatz [20], making it possible to derive the KPZ dynamical 
exponent z = 3/2 and various other quantities exactly. Other rigorous 
results concerning shape fluctuations in the ASEP can also be found in 
Ref. [21]. 

In order to formulate a mean field theory for the single step model, let 
N u (n, t) and N d (n, t) be the probabilities of finding an upward or downward 
segment with their lower edge rooted at height level n. Let us first consider 
a deposition process, in which a local minimum at level n is flipped into a 
local maximum at level n + 2. Having selected a random site, the probability 
to find such a local minimum at a given height can be approximated as fol- 
lows. Clearly, the probability of finding a downward segment on the left side 
terminating at height n is N d (n, t)/L, where L is the system size. With this 
probability, knowing that the height of the selected site is n, the adjacent seg- 
ment to the right can only go up or down so that the conditional probability 
to find an upward segment is given by N u (n, t) / (N u (n, t) + N d (n — 1, i)). Ig- 
noring possible correlations the total probability of finding a local minimum 
at height n is the product of these two expressions. The deposition process, 
taking place with probability p, therefore leads to a loss of probability at 
level n 

p N d {n,t)N u {n,t) 



N u (n, t) -► N u (n, t + dt) = N u (n, t) 
N d (n, t) N d (n, t + dt) = N d (n, t) 



L (N u (n,t) + N d (n-l,t)) 
p N d (n,t)N u (n,t) 
L (N u (n,t) + N d (n-l,t)) 



(9) 



3 



Initially the velocity is higher, and the excess velocity \v(t) — v^] decays as t 1//3 . 
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and a corresponding gain at level n + 1. Similar expressions can be derived 
for the evaporation process. Obviously, this approximation accounts for the 
restriction Q and the one-dimensional structure of the model but disregards 
possible nearest-neighbor correlations. 

The structure of Eq. © suggests that the probabilities N u (n,t) and 
Nd(n,t) evolve exactly in the same way. In fact, it is easy to see that in 
a system with periodic boundary conditions the numbers of upward and 
downward segments are exactly equal. Assuming the same to hold in an 
infinite system we have 

N u (n,t) = N d {n,t) (10) 

for every n and t. Thus, introducing a combined probability density 

N d {n,t) + N u (n,t) 



P(n,t) 



2L 



(11) 



the loss at level n due to a deposition event in Eq. © can be recast as 

Pnit? 



P(n, t) P(n, t + dt) = P(n, t) - p 



P n (t) + P„_i(i) 



(12) 



Collecting all loss and gain contributions due to deposition and evaporation 
one arrives at the following set of mean field equations 



P 



n-l 



Pn{tf 



+ (1-P) 



Pn-l(t) + Pn-2(t) 
Pn + l{t? 



P„(t)+P n _l(t) 
Pn{t? 



P n+l it) + P n+2 (t) P n (t) + Pn+1 (t) 



(13) 



which serve as a starting point for all further calculations throughout this 
paper. Notice that the form of the denominators appearing on the r.h.s of 
Eq. (|TH|) is a consequence of the restriction (jJJ), and that some care has to 
be taken when one of them vanishes. Since the numerators are quadratic we 
assume that each of these terms is zero whenever their denominator vanishes. 



Introducing a probability current flowing between neighboring levels 



J, 



n,n+l 



(t)=p- 



Pn(tf 



P n (t) + P n _i(t) 
these equations can be also written as 
d 



(1-p) 



P n+1 (t) + P n+2 {t) 



ot 



Pn{t) — Jn-l,n{t) ~ Jn,n+l(t) 



(14) 



(15) 
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Obviously, they conserve probability Y2k=-oo Pk (*) so that the integrated 
probability distribution 

n 

Q n (t):= Y, p k(t) (16) 

fc=— oo 

satisfies the simple evolution equation 

d 

g- t Qn(t) = -J n ,n+l(t). (17) 

By construction these mean field equations reflect both the one-dimensional 
structure as well as the restriction but they ignore spatial correlations be- 
tween the segments. The full model does exhibit such correlations, but it 
evolves towards a trivial state without correlations (corresponding to a sim- 
ple product state in the ASEP). Although this trivial state is never reached 
in an infinite system, it may explain why the mean field equations proposed 
here reproduce so many of the observed phenomena faithfully, some of them 
even exactly, as will be shown in the following sections. 



3 Exact solution of the mean field equations 



Let us first consider the case of a free interface, where the height index n 
runs over all integers from — oo to +oo. 

As the KPZ equation is invariant under appropriate rescaling of space, 
time and height [5], we can carry out the continuum limit by introducing a 
new height variable h = nA, where A is the new height unit of the rescaled 
system. In order to investigate the asymptotic properties of the roughening 
processes let us assume that P n {t) varies only slowly with n and expand the 
r.h.s. of Eq. (|13|) as a Taylor series around h. Keeping contributions up to 
fourth order in A we obtain the partial differential equation 



(18) 



A 3 (2p- 1) 
24 

A 4 {2P'(h,t 



3P'(h,tf 6P'(h,t)P"{h,t) 



P(h,t) 2 P(h,t) 
i 4 5P'{h,t) 2 P"(h,t) 



AP"'(h,t) 



+ 



P(h,t) 3 P(h,t) 2 
2P"(h,t) 2 2P'(h,t)P"'(h,t) 



+ 



P(M) 



+ 



P(h,t) 



P""(h,t) 
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where the prime stands for a partial derivative with respect to h. Obviously 
the leading term of order A on the r.h.s. generates a uniform propagation of 
the probability distribution, hence, the average height (h(t)) of the interface 
will asymptotically grow with the linear velocity 



plus some sub linear correction terms. Assuming ordinary Family- Vicsek 
scaling [3], it is therefore near at hand to test the validity of the scaling 
form 



Notice that the normalization of the height probability distribution implies 
the scaling function f(z) to be normalized as well. In what follows we solve 
Eq. ()18|) both in the equilibrium case p = 1/2 and the non-equilibrium case 
p ^ 1/2 confirming that our results do not depend on A. In particular, we 
will show that higher order terms appearing in the expansion of Eq. (|13|) 
turns to be irrelevant, vanishing in the asymptotic limit t — > oo. The correct 
asymptotic behavior of Eq. (jlMjl will be therefore recovered by setting A = 1. 

3.1 Equilibrium roughening of a free interface 

We start analyzing the special case p = 1/2, where the dynamic processes 
of the full model are known to exhibit detailed balance. In this case the 
velocity v is zero and the first-order and third-order contributions on the 
r.h.s. of Eq. ()18|) vanish. Inserting the Ansatz ()20|) into Eq. (|18|) we find 
that, up to fourth order, the partial differential equation reduces to a non- 
trivial ordinary differential equation for the scaling function if and only if 
7eq = 1/4 (the subscript denoting the equilibrium case) . This is exactly the 
value predicted by the EW theory for equilibrium roughening. Moreover one 
easily notices that by fixing 7 eq = 1/4 higher order terms 0(A 5 ) occurring 
in the Taylor expansion of Eq. (|13j) are irrelevant in the asymptotic limit 
t — ► oo. The differential equation therefore reads 




(19) 




/eq(^) 3 



1 



/ eq (z) 4 - ^f eq (z)f^(zff^(z) + 

A 4 / eq w 2 (Aw 2 +/;w/:^)) + 



(21) 
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where z = ht 7cq denotes the scaling variable. Integrating both sides we 
obtain 



A 4 



^feq( Z )fcq( Z J , ,,/// 



feq(z) 



This equation admits the two simple solutions 



f:: e {z) 



z 



1 ' --' 



: CXp 



and 



^bound 
' eq 



(*) 



2 l/4 



z 2 exp 



■— 1 



(22) 



(23) 



(24) 



which have been normalized over the real line. The first solution /gq 6e is 



a simple Gaussian and represents the physical solution for a free interface 
starting with a flat initial condition h(x, t) = 0. The second solution /,J^' und 
is characterized by two different maxima over the real line and is therefore 
dismissed as unphysical in the free case. However, as we will see in Sect. 4, 
this solution becomes physically meaningful in the presence of a hard-core 
wall at zero height. 

To summarize we note that the mean field equation for p = 1/2 does 
indeed capture the features of one-dimensional EW roughening in the ther- 
modynamic limit L — » oo. 



3.2 Nonequilibrium roughening of a free interface 

We now turn to the nonequilibrium case p ^ 1/2. Inserting again the scaling 
form ()2U|) and the expression for the velocity (|19|) into the partial differential 
equation (|18|). we find that, up to third order, the partial differential equation 
reduces to a non-trivial ordinary differential equation for the scaling function 
(i.e., without explicit occurance of t) if and only if 7 = 1/3. The differential 
equation then reads 

' \8f(zf + 3kf'(zf - 6kf(z)f'(z)f"(z) + (25) 



f{zf 



Af(zf{2zf\z) + kf"{z)) 



0. 



where z = (h — vt) t 7 denotes the scaling variable in the comoving reference 
frame and 

k = (2p- 1)A 3 ^ 0. (26) 
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As in the equilibrium case the postulate of Family- Vicsek scaling applied 
to the mean field equation already determines the roughening exponent. 
Remarkably the value 7 = 1/3 coincides exactly with the known value of a 
KPZ process in 1+1 dimensions. 

As can be verified easily, upon fixing 7 = 1/3, the fourth order terms 
(and all higher order terms) of the Taylor expansion turn out to be irrelevant 
in the asymptotic limit t — > 00 and thus generate only short time corrections 
to the scaling function. It is also worth commenting that asymptotic EW 
scaling can only be seen in the symmetric case p = 1/2. For any small devia- 
tion from this value the third-order terms do not vanish, leading eventually 
to a crossover to KPZ scaling in the limit t — > 00. Therefore, the mean 
field equations nicely reproduce the character of the KPZ nonlinearity as a 
relevant perturbation. 

Assuming that f(z) 7^ and integrating both sides of Eq. (|25|) one 
obtains a simplified equation which, by substituting f(z) = u(z) 4 , can be 
further reduced to a simple Airy differential equation 



where Ai(z) and Bi(z) are Airy functions (see for instance [22]). For given 
\k\ the two solutions differ only by a reflection z — > — z so that for the rest of 
this section we can restrict ourselves to the case k > (which corresponds 
to a negative nonlinear term, i.e. A < 0). 

The two integration constants have to be chosen such that f{z) is prop- 
erly normalized and the appropriate boundary conditions are satisfied. For 
a free interface, the scaling function f(z) has to vanish for z — * ±00 in such 
a way that all of its moments are finite, hence C2 = 0. Surprisingly, the 
remaining solution oscillates for z — > 00 and does not vanish fast enough 
to yield finite moments. We conclude that the physically meaningful solu- 
tion extends from z = — 00 to the first root of the Airy function zq (where 
f(zo) = f'(zo) = f"(zo) = 0) and vanishes elsewhere. The solution for the 



zu(z) + 2ku"{z) = 



(27) 



with the general solution 




for k > 



for k < 



(28) 
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Figure 4: Rescaled skewed profile obtained within the mean field approximation 
(solid line) compared to the numerically determined KPZ profile using the single 
step model with p = 1. Differences between the two profiles can be highlighted 
by their mean value a\. The mean value controls the KPZ excess velocity and it 
is known to scale as A 1 / 3 . While the mean field profile is characterized by a\ ~ 
1.13184, direct numerical estimate renders a\ t ssM — 0.60. The vertical axis is 
plotted in a logarithmic scale. 



free interface therefore reads 

Ai 





/(*) 



i 

AT 



(2k) 1 / 3 



for — oo < z < zq 

for Zq < Z < +00 



(29) 



where z ~ 2.94583 /c 1 / 3 and M ~ 0.127153 A; 1 / 3 . As it can be seen, the 
parameter k, describing the strength of the KPZ nonlinearity, appears here 
as a simple metric factor in the scaling function. Notice that by Eq. (|26|) 
the height unit A has been absorbed in k. 

Figure 0] shows the solution (|29|) in comparison with the numerically 
determined profile of a freely roughening KPZ interface. Although the two 
curves are different due to the approximative character of the mean field 
theory, they share essential qualitative features. To quantify those it is 
instructive to compute the skewness 



S 



_£3_ 

3/2 

c 2 



03 — 3a2ai + 2a 3 



(02 



l?)3/2 



(30) 



where c n is the n th moment of the height probability distribution and a n = 
f*° dz z n f(z) denotes the n th central moment of the rescaled profile f(z). 
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The skewness is expected to be universal in modulo for all KPZ growth 
processes, with its sign being equal to the sign of the nonlinear term. Known 
numerical estimates [23] give the value |Skpz| = 0.28±0.04, which is in good 
agreement with the theoretically computed skewness for the PNG model 4 [8], 

Spng - 0.2935 . (31) 
Mean field theory, on the other hand, renders the value 

Smf ±0.465970 , (32) 

where the positive (negative) sign correspond to the case k < (k > 0). 
Although this value is different from direct numerical estimates, it has the 
correct sign and the same order of magnitude, showing that the mean field 
theory captures qualitatively the influence of the nonlinear KPZ term. 

Surprisingly, the mean field theory predicts that the interface profile is 
asymptotically bounded for z > at a finite value zq. This means that within 
mean field the advancing front of the distribution exhibits a sharp cutoff 
rather than a stretched exponential tail. However, on the opposite side, 
where z is negative, the profile does indeed decay as a stretched exponential: 

f{z) ~ \z\-V 4 e~^ ]z]3/2 (z-»-oo) (33) 

This result suggests that r] + = 3/2, which coincides with the theoretical 
value predicted in the context of directed polymers and of the PNG model. 
On the other hand, r/_ does not exist within the mean field approximation, 
which therefore fails to correctly describe large negative (w.r.t. the sign of 
the KPZ nonlinearity) height fluctuations. 

4 Roughening in the vicinity of a wall 

We now modify the single step model and the associated mean field equa- 
tions in order to incorporate the presence of a hard core wall. Our aim is to 
determine the surface exponent 6 introduced in Eq. @ within the mean field 
approximation. In terms of the continuous height variable h, the density of 
pinned sites can be defined as the integral of the height probability distri- 
bution between the hard core wall and some arbitrary small height level ho, 

4 The PNG model is characterized by a positive nonlinear term. 
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i.e., 

rha+vt rh t T 

p(t) = / P{h,t)dh= / , (34) 

Jvt JO 

so that the surface exponent is completely determined by the behavior of 
the scaling function f(z) for z -C 1. 

In order to formulate the appropriate boundary condition in the con- 
tinuum limit, one has to resort to the discrete formulation of the problem. 
Following the approach outlined in Ref. [24], the wall is initially located at 
zero height and moves dis continuously with the average velocity v = p — | , 
its actual height being given by n (t) = [vt\ (where |_-J indicates the integer 
part). The interface is restricted to evolve above the wall, i.e. 

Pn{t) = (n < no(t)) (35) 

so that the mean field equations ()13|) have to be modified accordingly close 
to the wall. In particular there is no probability current between level no 
and no + 1, so that 

d 

Q^Pno+lyt) = -Jno+l,no+2{t) , (36) 

while Eq. qi5|) still holds for n > no + 1. Depending on p one has to distin- 
guish three different cases. If v > the wall advances by one unit in time 
intervals At = 1/v , flipping up all local minima at level no + 1 by two units. 
This means that P no +i is increased by P no while P no is set zero. If v < the 
wall retracts by one unit in time intervals At = \l/v\, allowing height level 
no, which was previously set to zero, to become nonzero during the subse- 
quent evolution. Finally, for v = the wall does not move, i.e. no = for all 
times t. The moving wall makes it difficult to specify the correct boundary 
conditions, so out of equilibrium we will derive them in the special cases 
p = 1 and p = 0, where the KPZ nonlinearity is maximal. Our reasoning, 
which once more relies on a series expansion in the proximities of the wall, 
shows that both a pushing (p = 1) and a retracting [p = I wall impose a 
Dirichlet boundary condition for the scaling function. Surprisingly, it turns 
out that a retracting wall (p = 0) does not fix any boundary condition for the 
scaling function, which is free to assume any finite value at wall level, thus 
justifying the high density of pinned sites which is numerically observed in 
the case A > (see Fig. |2j). General scaling arguments suggest that results 
obtained for p = 1 (p = 0) hold for any p > 1/2 (p < 1/2). 



16 



Figure 5: Rescaled mean field profiles in the presence of a hard core wall compared 
to the corresponding profiles in the SSM for p = 1 (left panel) and p — (right 
panel). The vertical axis have been plotted in a logarithmic scale. 

4.1 Depinning in the equilibrium case p — 1/2 

In the equilibrium case we have v = so that the wall does not move. As 
shown in the Appendix, a wall at zero height imposes a Dirichlet boundary 
condition /(0) = 0. Obviously, the only solution satisfying this boundary 
condition is Eq. (|24|) 



which has been normalized here over the positive real axis. With this so- 
lution we can immediately read off the surface exponent from Eq. (|34|) . 



We note that this value coincides exactly with the known exponent for Ed- 
wards Wilkinson growth next to a wall. 

4.2 Depinning in the non-equilibrium case p = 1 

For p > 1/2 the wall advances discontinuously which makes it more difficult 
to specify the boundary conditions. As shown in the Appendix, in this case 
the co- moving wall again leads to a Dirichlet boundary condition /(0) = 
for the scaling function. According to Eq. (|28j) the corresponding solution 




(37) 



« = 3/4. 



(38) 
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then reads 



/(*) 



i 

77 




Ai 



(2k) 



1/3 



^Bi 



(2k) 1 / 3 



for < z < Z\ 
for zi < z < +oo 



(39) 

where z\ ~ 3.32426 A; 1 / 3 is the first positive root of the scaling function f(z) 
(where also its first three derivatives vanish) and M ~ 0.133454 A: 1 / 3 . As in 
the free case the profile exhibits a sharp cutoff (see Fig. [5J), although at a 



different value of z. Since f{z) 
by 



~ z 



for z — > 0, the surface exponent is given 



oMF 



4/3, 



This values has to be compared with the numerical estimate 
in Eq. ©. 



(40) 
1.1(1) 



4.3 Depinning in the non-equilibrium case p = 



For p < 1/2 the wall moves discontinuously backward. As shown in the 
Appendix, this situation is special in so far as the retracting wall does not 
specify any boundary condition on the scaling function, allowing /(0) to be 
nonzero. In fact, according to Eq. 1)28(1 the only normalizable solution is 
given by 

'< 2 > " w 



Hw^>\ (41) 

with the normalization constant M ~ 0.00584355 A; 1 / 3 . Since /(0) > the 
surface exponent is given by 



oMF 
7 p=0 



1/3, 



This values has to be compared with the numerical estimate 
in Eq. ©. 



(42) 
0.22(2) 



5 Conclusions 

In this paper we presented a mean field theory for 1+1 dimensional nonlinear 
growth processes evolving according to a KPZ equation. It is worth stressing 
that our approach does not neglect all types of fluctuations, as the term mean 
field usually suggests. Instead our equations retain information about the 
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one dimensional spatial structure and the height restriction hi — hi+i = ±1. 
Therefore, our mean field theory is not expected to hold exactly above some 
upper critical dimension, rather it serves as an approximation for the one 
dimensional case. Although we neglect spatial correlation between local 
slopes of the roughening interface, the theory has a very predictive power. 
Its success can be ascribed to a fluctuation dissipation theorem which is 
known to hold only for the 1+1 dimensional case [5]. Moreover, the mean 
field equations presented herein have been derived in the thermodynamic 
limit L — > oo, and thus they are suited for probing finite time behavior, i.e. 
t<^L z . 

Our approach is therefore successful in predicting a power law decay for 
the density of interfacial sites pinned to the substrate, thus supporting previ- 
ous numerical studies of the non-equilibrium case. Although the mean field 
surface exponents 6 differ from the numerically estimated values, our theory 
correctly reproduces the dramatic difference between the surface behavior 
in the presence of a negative or positive nonlinear KPZ term. 

For what concerns the finite-time bulk properties of a KPZ interface, 
our simple mean field theory correctly reproduces the exact value for the 
roughening exponent 7 = 1/3 and the skewed nature of the finite time 
height probability distribution. We found that one of the two tails of such 
a distribution decays as a stretched exponential with exponent r] + = 3/2, 
thus confirming previous results obtained in the context of directed polymers 
and for the PNG model. While our theory successfully predicts large positive 
(w.r.t. to the sign of the KPZ nonlinearity) height fluctuations, it fails in 
describing large negative ones, exhibiting a sharp cut-off for the negative tail. 
It is therefore interesting to note that a simple scaling argument [9] proposed 
in the context of the PNG model directly relates 77+ to the roughening 
exponent, i.e. rj + = 1/(1 — 7), while no corresponding argument can be 
worked out for ry_. This is due to the fact that height fluctuations with 
the same sign w.r.t. the KPZ nonlinear term manifest as "bumps" (or 
"holes") which grows laterally, while height fluctuations with the opposite 
sign manifest as "holes" (or" bumps" ) which shrink laterally. 

It is also worth noticing that the mean field theory reproduces correctly 
almost all bulk nonlinear critical properties at finite times, while it gives 
only approximate results for the surface exponent 8. This is an indication 
that the substrate introduces spatial correlations between local slopes at low 
height levels. It is our belief that a detailed study of the SSM with a hard 
substrate may eventually lead to the exact analytical knowledge of critical 
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depinning properties. It would also be interesting to find out whether the 
methods introduced in Ref. [8] can be applied to the problem of a KPZ 
interface with a wall. 

Finally, equilibrium results are correctly reproduced as a marginal case, 
and small out of equilibrium corrections eventually leads to full KPZ behav- 
ior. 
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A Boundary conditions in the presence of a hard 
core wall 



In order to solve the mean field equation Q18|) when a hard core wall is 
imposed, it is necessary to go back to the discrete formulation of the problem 
and to consider the proper evolution equation for the height probability 
distribution ()36|) close to the wall. In the following we analyze the special 
cases v = and v = ±| in order to derive the corresponding boundary 
conditions for the scaling function f{z). 



A.l p = |, v = 

In this case the wall does not move and no = for all times t. According to 
Eq. 1)36|) the density of contact points P\(t) then evolves as 

^Pa(i) = -J lj2 (i) , (43) 

where 

(f) _Pi(t) P 2 (t) 2 

In the limit t — > oo this equation implies a Dirichlet boundary condition 
/(0) = 0. To see this let us assume that /(0) ^ with |/'(0)| < oo. 
Assuming EW scaling the first three probabilities would be given by P\(t) — 
P2(t) — Ps(t) — i -1//4 /(0) (where we expanded the scaling function / around 
z = keeping only the leading term), giving rise to a current Jizit) — 
i^-i/ 4 y(0). Since the l.h.s. scales as i _5//4 , Eq. (|4^)) cannot hold unless 
/(0) = 0. 
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A. 2 p=l,v = l 



For p = 1 one has At = 2, i.e. the wall advances by one unit after every 
second time step. To find the appropriate boundary condition for f(z) (now 
assuming KPZ scaling with 7 = 1/3), let us first consider the continuous 
temporal evolution between two advancements when the wall is fixed at 
some height no[— to/2\ . As in the previous case the probabilities P n (t) with 
n < no vanish. According to Eqs. (|14jl. (|T5|) and (JBEJ), the height probability 
distribution at the first two levels above the wall obeys to the differential 
equations 

d 

wrPn +i{t) = --Pno+H*) 

d > p2 Q+2 (45) 

^7-Pn +2(i) = +P no +l(t) ~ — • 

Let us again suppose that /(0) 7^ and |/'(0)| < 00, i.e., just after the ad- 
vancement of the wall at time to we assume that to leading order P no+ i(to) ~ 
Pn +2{to) — t Q /(0) =: c(to). Iterating the differential equations (|loT) over 
two time steps, and by assuming c(t) ~ c(to) for t$ < t < to + 2, one obtains 
Pn +i(to + 2) — c{to)e~ 2 . On the other hand, by numerically solving the 
differential equation for level uq + 2 

P no+2 (t) ~ +C (t )e-^ - *?+ 2 (46) 

c(t )e (.* *oj + P no+2 

one gets P no+ 2(to + 2) « 0.596 c(to)- At time to + 2 the wall advances by 
one unit, i.e. all local minima at height no are flipped upwards, meaning 
that P no+ i(to + 2) is first added to P m+ 2(to + 2) and then set to zero. Just 
after advancement P rao +2(to + 2) ~ 0.732 c(io) which is in contradiction the 
assumption unless /(0) = 0. Hence for p = 1 the wall imposes again a 
Dirichlet boundary condition. 



A.3 p = 0,v = -± 

For p < the wall retracts by one unit after every second time step. Let us 
first consider the continuous temporal evolution between two moves when 
the wall is fixed at height no = [—to/2\. As usual the probabilities P n (t) 
with n < no vanish and the he first two levels, where the height probability 
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distribution is nonzero, evolve according to the differential equations 



fi P 2 

° p (+\ _ r n +2 



a p 2 p 2 1 ' 



9t Pn +2 + Pn +3 Pn +3 + Pn +4 

Just after retraction level uq + 1 can be visited by the interface for the first 
time, hence P no +i(to) is initially zero and becomes nonzero as time evolves. 
Assuming that /(0) 7^ and |/'(0)| < 00, to leading order the other prob- 

— 1/3 

abilities have the initial values P no+ 2(to) ~ P no+ ^(to) ~ t /(0) =: c(to)- 
Iterating the differential equations (|45jl over two time steps one obtains 
P no+ i(to + 2) ~ c(to) while the higher levels remain unchanged. Thus the 
almost constant probability distribution in the vicinity the wall is simply 
'extended' to the new level that becomes available by retraction of the wall, 
meaning that a nonzero value /(0) 7^ and |/'(0)| < 00 is consistent with 
the equations (|47j). Loosely speaking the wall retracts so quickly that it does 
not impose a specific boundary condition, allowing /(0) to take any positive 
value. The equations are built in such a way that this value is simply copied 
from the following height level. 
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